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Abstract. We present a detailed numerical study of the dynamics of a disordered one- dimensional Bose- 
Einstein condensates in position and momentum space. We particularly focus on the region where non- 
linearity and disorder simultaneously effect the time propagation of the condensate as well as the possible 
interference between various parts of the matter wave. We report oscillation between spatially extended and 
localized behavior for the propagating condensate which dies down with increasing non-linearity. We also 
report intriguing behavior of the phase fluctuation and the coherence properties of the matter wave. We 
also briefly compare these behavior with that of a two-dimensional condensate. We mention the relevance 
of our results to the related experiments on Anderson localization and indicate the possibility of future 
experiments 

PACS. 03.75 -b Matter waves - 37.10. Jk Atoms in optical lattice - 05.60.-k Transport processes 



1 Introduction 

The spatial behavior of the time dependent density profile 
of a cold atomic Bose-Einstein condensate (BEG) in pres- 
ence of a optical disorder potential has been under intense 
experimental [1], [2], [3], [4] and theoretical investigation 
[S], [5], [3, [Sj, [in]in recent years. Such interest in this 
issue is attributed to the possibility of observing the phe- 
nomenon called Anderson localization (AL) ^llj of matter 
waves. AL corresponds to the spatial localization of a wave 
due to destructive interference between multiply scattered 
waves from the modulations of a disorder potential. The 
BEG of ultra cold atoms is considered as an ideal system 
to observe such localization. This is because of the system 
forms a giant matter wave of ultra cold bosonic atoms with 
very weak interaction among them. Additionally the sys- 
tem can be prepared in reduced dimension and according 
to the scaling theory of localization [T^] the effect of disor- 
der gets enhanced in that case. In two recent experiments 
carried out in Orsay [13] and in LENS [T^ Anderson lo- 
calizations of matter wave were observed respectively in a 
weakly interacting quasi one dimensional repulsive con- 
densate and non-interacting one dimensional condensate. 
The rapid progress in this field also has been recently sum- 
marized in two recent reviews by Orsay and LENS group 

This paper investigates the the joint effect of non- 
linearity and disorder on the time evolution of the co- 
ordinate and momentum space density, phase fluctuation 
and coherence properties of such condensate We begin 
with the study of time dependent density distribution in 
co-ordinate as well as in momentum space. According to 



the scaling theory of localization 12J in one spatial dimen- 
sion the motion of a quantum mechanical particle is either 
ballistic, or localized in presence of disorder with infinites- 
imal strength, and the diffusion is absent. However recent 
analysis [8], [9], [15] showed that the presence of interaction 
changes this situation and localization properties become 
momentum dependent. Different wave- vector components 
of the propagating Bose-Einstein condensates will behave 
differently as they get scattered from the random modula- 
tions of the disordered potentials. And this motivates us to 
numerically investigate the localization properties of the 
fast as well as slow components of the static condensate. 
Additionally we study the phase fluctuation and the asso- 
ciated coherence of the condensate as a function of time 
which provides new and interesting informations about 
this system. In lower dimensional systems strong phase 
fluctuations quite often destroys the long range phase co- 
herence in Bose-Einstein systems. Though it was shown 
[17] that for a low dimensional condensate the phase coher- 
ence exists over a length scale comparable with the typical 
system size, for the time evolved condensate in presence 
of a disordered potential this argument cannot be auto- 
matically extended. This motivates us to investigate the 
coherence properties of such time evolved condensate at 
least numerically. 

To place our analysis in the experimental context let 
us mention that the issue addressed in most of the cur- 
rent experiments and the associated theoretical studies is 
to ascertain under what condition the system will mani- 
fest AL and what will be the impact of non-linearity on 
this localization process. The AL has been confirmed by 
the presence of an exponentially decaying envelope wave- 
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function of the relevant system. In a recent experiment [13 
the time of flight image of the edge of a quasi-one dimen- 
sional Rubidium BEC which has low density and almost 
negligible interaction energy, showed this exponential de- 
cay. Thus it has been concluded as Anderson localized. 
In this experiment both the strength of the disorder as 
well as the density of the condensate are kept at a rel- 
atively lower value to achieve the regime of AL. These 
experimental conditions are borne out of the fact that a 
very strong disorder classically localizes a wave purely on 
energetic ground whereas for a high density condensate 
cloud the interaction effect is too strong to allow the An- 
derson localization to happen. The preceding theoretical 
analysis [8], [9] which has identified these experimental 
conditions suitable for the observation of Anderson local- 
ization can be summarized as follows. A typical trapped 
condensate contains a range of momentum components. 
The higher momentum components that mostly reside at 
the low density edge of the condensate moves much faster 
as the condensate propagates through the disorder poten- 
tial. After an initial stage of expansion in which most of 
the interaction energy of the condensate gets converted 
into kinetic energy, these higher momentum components 
of the condensate are mostly present at the edge of the 
expanded condensate and behave like plane wave. These 
are almost free from the non-linear effect and thus are 
appropriate for observing AL. 

This should be contrasted to the behavior of the cen- 
tral part of the time evolved condensate which is mostly 
populated by the low momentum components and con- 
tains non-linear effect. The time evolved density profile of 
this central region nevertheless gets strongly modified by 
the disordered potential [15]. However, why such modifi- 
cations do not lead to AL is still not completely under- 
stood. On the other hand it is the behavior of this section 
of the cloud which manifests the interplay between the 
non-linearity and the disorder and sheds light on the im- 
pact of the non-linearity on the AL. Our work address this 
unresolved issue by trying to find out the relative impact 
of non-linearity and disorder by studying numerically the 
time evolution of density modulation, phase modulation 
and coherence properties which combines the effect of den- 
sity and phase. Thus this analysis provides us important 
information about the fate of interference induced quan- 
tum localization in presence of non-linearity and augment 
a number of theoretical studies that studied these issues 
previously [6] , [8] , [9j , [10] . 



[I3] 2.1 Model of the disorder potential 



The results that we shall show in the subsequent sections 
have been generated with the disorder potential as de- 
scribed in the work of Huntley [19] and has been used 
in a number of experimental work to model the disor- 
der potential. Here, the real and imaginary part of com- 
plex amplitude of light are constructed by Gaussian ran- 
dom numbers. Such a disorder potential is known as the 
speckle potential and the height of the speckle which char- 
acterizes the strength of the disorder potential is propor- 
tional to the standard deviation of the disorder potential 
around an almost zero mean value. We have compared the 
results obtained in presence of such disordered potential 
with ones where we have created the disordered potentials 
simply from a uniformly distributed random numbers [lOj 
or real random numbers with Gaussian distribution. The 
later type of random potentials have non-zero mean and 
its strength is characterized by the first as well as the sec- 
ond moment of the distribution of the random variables. 
The details of the time dependent as well as stationary 
density profiles of condensates indeed differ from one real- 
ization of the disorder potential to another, but the basic 
features related to the localizations remained unaltered. 
All the results presented here are for single realization of 
the disorder potential. But repeating the calculation for 
various realizations of the disorder potential we found that 
the basic conclusions, these results deliver, do not depend 
on the details of a realization. 

Also, the mean field description of the disordered con- 
densate through the disordered Gross-Pitaevskii equation 
is valid ( for a more detailed explanation, see [15] ) as long 
as the disorder potential fluctuates over a length scale that 
is larger than the so called healing length of the conden- 
sate. To ensure that this condition is obeyed we discard 
from the Fourier spectrum of the disorder potential all 
wavenumbers that are above a given cutoff fee ~ 27r/Ac 
[To] .The inverse Fourier transform provides a random po- 
tential that varies on length scales larger or equal to Ac 
and which can be formally written as 
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where Vg is the original speckle potential. In the final cal- 
culation we have used this smoothened disorder potential. 



2.2 Mean field description of the disordered 
condensate 



2 The Theoretical framework 



In this section we discuss briefly various details of the 
theoretical framework of the problem. First we talk about 
the model disorder potential. Then we shall discuss the 
disordered non-linear Schrodinger equation which will be 
solved to study the behavior of the disordered condensate 



We have described the time evolution of the quasi-one di- 
mensional repulsive condensate through the mean field 
Gross-Pitaevskii (GP) equation. Most of the recent ex- 
periments [T]l3ll4lfT3] can be well explained by such mean 
field theory. From the pure theoretical perspective this 
also helps to address the issue if one dimensional Nonlin- 
ear Schrodinger equation (NLSE) shows Anderson local- 
ization in presence of a disorder potential [2U], that its 
linear counterpart does. 
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Mean field GP equation is obtained from the micro- 
scopic description of the ultracold bosonic systems by re- 
placing the bosonic field operator by its vacuum expecta- 
tion values which is taken as condensate order parameter. 
The details are given in a large number of text books ( for 
example, see ref. [21]) Typically for these ultra cold atoms 
the interaction is extremely short ranged and character- 
ized by a single parameter which is the s-wave scattering 
length. In the GP equation the interaction is characterized 
by a nonlinear term that is cubic in the order parameter. 

Our parameters are chosen such that we are close to 
the typical experimental regime [13j . A quasi one dimen- 
sional condensate is achieved by tightly trapping the con- 
densate in transverse direction such that the only quan- 
tum mechanical state occupied in the transverse direction 
is the ground state of a two-dimensional harmonic oscil- 
lator. The mean field Gross Pitaevskii equation that de- 
scribes such one dimensional condensate in presence of a 
speckle potential Vs is given by [T7] 

where gid = 2at)w± is the one dimensional effective cou- 
pling constant [18' that depends on the s-wave scatter- 
ing length a and the radial harmonic trapping frequency 
oj^ of the condensate. Here uiz is the trapping frequency 
in the longitudinal direction. We first rewrite the Gross- 
Pitaevskii equation in a dimensionless form in the same 
way as in [lOj . 
It becomes 

.di)^ 1 d'^ip^ 1 2 , , T. , , An , |2 , to\ 
= + 2^ Vw + Vdi^w + 9N\%pyj\ ipw (3) 



by choosing hw^ the dimension of energy, y j^j- as the 

unit of length, lo^^ as the unit of time. Vd is defined in 
Eq. [1] Since we shall be considering the time evolution of 
the condensate density after its release from the trap we 
set the trapping potential to zero. We also denote the 
dimensionless co-ordinate , ^ as x from now on where 

as for the time we set ujzt = r. ip^ that gives the conden- 
sate wave function in presence of disorder is normalized to 
unity in the above equation and the interaction coupling 
is thus proportional to the number of particles N . As dis- 
cussed in the section 12. 1[ the validity of this mean field 
equation in the presence of the disorder potential requires 
that the Ac > ^, where ^ is the healing length of the con- 
densate ( \J 2^ from Eq. [3]) . In the following numerical 

calculations S,highest = 0.250 for dimensionless chemical 
potential /i = 6, where as the Ac = 1.25, and, satisfies the 
above criterion. 

Using the parameters for ^"^Rb we found that under 
typical experimental condition 17 « .03 [13j . The values 
chemical potentials fi (in the unit of hcoz) that have been 
used in the numerical analysis pertains to ~ 10"^ — 10^ 
number of particles for a condensate. In this regime the 
Thomas- Fermi (TF) approximation works well [TTlfTU] and 
this is also close to the typical experimental conditions. 



We assume that at t = the condensate has a TF pro- 
file for a given chemical potential /i [10]. We have defined 
Ltf = -s/^/I as half of the size of such TF condensate. 
Thus the total number of particles N, that can be obtained 
by integrating TF density over this TF length, increases 
with increasing /i. Since TV appears as the coefficient of 
the non-linear term in eq. [3l by increasing /i one actually 
increases the strength of the non-linearity and make the 
condensate more interacting. 

We shall obtain the condensate wave function and the 
corresponding density, pw{x,t) = |-0iu(a;, at a later 
time t by numerically integrating the above time depen- 
dent equation ([3]) in the method outlined in [22] and has 
been used in jlOj over a discrete grid whose linear size is 
much larger than the TF length of the condensate. We 
have also tried our numerical code for several grid sizes 
to make sure that the discretization process does not add 
any spurious feature to our result. We will discuss more on 
the effect of numerical discretization in the next section 
with the help of results. Additionally we have also used an 
absorbing boundary condition so that any possibility of re- 
flection from the boundary is eliminated. This describes 
the theoretical framework and the values of parameters 
that are used in our numerical analysis. In the following 
section we shall present our results and will analyze them. 

3 Results and discussions for one dimensional 
condensate 

3.1 Time dependent density 

It has been shown [23] that for a propagating condensate, 
most of the interaction energy gets converted into kinetic 
energy over the time t^-^^Ldzt = T^\. Here we 
present the result at r = 10 in the Fig.[l]and|2] Each cal- 
culation have been performed for a number of realization 
of the disorder potential and features common to all such 
realizations are demonstarted graphically using one such 
representative realization. From time to time, we shall also 
point out those features which may vary from one realiza- 
tion of the disorder potential to another. 

All the subsequent figures broadly answer the following 
three questions : 

1. How the density, phase and the coherence properties 
of the condensate evolves when the disorder potential 
is kept same but the strength of non-linearity (/x)is 
varied ? 

2. How the above properties evolve with time when the 
non-linearity is kept same, but the strength of the dis- 
order potential is varied ? 

3. Finally, how the time evolution of all these properties 
changes when the higher Fourier components are re- 
moved from the initial profile of the condensate?This 
removal is done by only allowing the components sat- 
isfying ULtf < 1 at i = in the profile of the con- 
densate where Ltf is the TF length of the static con- 
densate. Henceforth, to avoid repetition, in the rest of 
the discussion we shall refer this wave-number filtered 
condensate profile as slow condensate. 
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Fig. 1. color online: one dimensional condensate in a disor- 
der potential. The disorder potential has mean and standard 
deviation 0.59. Ac for the disorder potential is 1.25. All quanti- 
ties are in the dimensionless unit, (a) The left column gives co- 
ordinate space density profile density profile of an one dimen- 
sional condensate after time of flight of uj^t — 10 for ^ — 10 
and /X = 15. x axis is in unit of TF length Ltf of the static 
condensate in each case. The number of particles varies from 
10^ - lO'' as the chemical potential changes (b) Corresponding 
time evolved density profile of the "slow condensate such that 
at t = 0, kLxF < 1 are given in the right column for the same 



Whereas the importance of first two questions can be eas- 
ily appreciated the relevance of the last question demands 
some discussion. In reference [9| for the case of a two- 
dimensional condensate with radially symmetric density 
profile it has been argued that that the components which 
satisfies the criterion kLxF < 1 will show a diffusive a be- 
havior after the initial stage of ballistic expansion. In that 
case k is the absolute magnitude of the two dimensional 
wave vector and Ltf corresponds to TF radius. In one di- 
mension, in the non-interacting limit no diffusive motion 
can take place. However, it can be generally expected that 
the faster (kLxp > 1) and slower component will response 
to the disorder diflFerently. Reference [15] also describes in 
detail how this leads to the existence of a mobility edge 
in one dimension in presence of non-linearity. It is in this 
context the above mentioned numerical investigation has 
been carried out. We find that slowing down of the con- 
densate in this way indeed effects its transport through 
the disorder. 
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Fig. 2. color onhne: Momentum space profile of one dimen- 
sional condensate in a disorder potential. The disorder poten- 
tial as well as other parameters are same as in the case of Fig. 
\j](a) The left column provides the momentum (k) space profile 
of the actual condensate, (b) The right column provide the k 
space density profile for the "slow condensate" 



The column (a) of Fig. ([T]) shows that density profile 
for the full condensate as the chemical potential is being 
increased for the same realization of the disorder potential. 
Unless otherwise stated in all the figures chemical poten- 
tial is in the unit of fkOz- As pointed out, higher /i implies 
higher density of particles and thus higher strength of non- 
linear term. The left column shows a strong localization 
as the non-linearity is increased by changing fj, — 10 to 
fj, = 15. The lower plot in column (6) shows for the corre- 
sponding slow condensate shows that the removal of the 
higher k density components from the initial profile resists 
this strong localization at /i = 15. The patches that ap- 
pear in the plots of the condensate density in Fig. [T] are 
actually very fine density modulations which can be seen 
by enlarging these portions of these figures. The edge pro- 
file of all these disordered condensate shows a decaying 
envelope, where as bulk density modulation does not have 
such decaying envelope for the bulk portion. 

Given the fact that we are considering a repulsive con- 
densate, this strong localization for higher non-linearity is 
a surprising result. However it may be noted that in the 
context of discrete non-linear Schrodinger equation which 
is used to model a dilute BEG in presence of an additional 
periodic potential under the tight binding approximation, 
localized excitations due to both discreteness as well as re- 
pulsive non-linearity [24] was also reported. Thus one may 
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argue whether this locahzation is an artifact of the numer- 
ical discretization of the otherwise continuous functions in 
the GP equation. We will further discuss the possible ef- 
fect of numerical discretization with the help of fc-space 
density distribution in the next paragraph. 

To complement the above finding, in Fig. [2] we plot 
the fe-space density profile of the condensate for all the 
situations depicted in Fig.[T] For fi — 10, the momentum 
space density distribution of the time evolved condensate 
remains strongly localized around fc = with two very 
small peaks at the corner of the first Brillouin zone. As 
it is increased to /i = 15, we see that density peaks occur 
at fc = as well at the two corners of the first Brillouin 
zone. We note that we have considered in the numerical 
method the fc-space distribution within the first Brillouin 
zone only. In the right side of the same figure, namely col- 
umn (b) , we see that the removal of the faster components 
from the initial condensate density leads to a localization 
of the time evolved condensate density at different points 
along the fc-axis which are far removed from fc = and are 
located near the first Brillouin zone corners. The scatter- 
ing at the Brliiouin zone corners are a signature of Bragg 
scattering since kgcatter — kinuiai = G where G is the re- 
ciprocal lattice vector of the lattice where the scattering 
takes place. Here the only source of the periodic lattice 
is the numerical discretization. However numerical grid 
remains same throughout our calcualtion whereas the fc- 
space density Pw{k) is not the same at the Brillouin zone 
corner for various values of fj, as well as for the various 
strength of the disorder potential. It also changes sub- 
stantially for the slow condensate ( Fig. [21 column (b) ) 
Moreover we checked that for other size of the numerical 
grid various features of the result remains almost same. 
Thus we may conclude that the density peaks at the cor- 
ner of Brillouin zone is partly due to the the numerical dis- 
creteness and the artificial lattice which leads to a strong 
Bragg scattering. Additionally, the peak also depends on 
the strength of non-linearity and disorder and also on the 
initial momentum distribution. However within the frame- 
work of this numerical investigation we are unable to find 
out to what extent the Bragg scattering is influenced by 
the later set of effects. 

In Fig. [3] we have compared the above mentioned den- 
sity profiles with the case when there is no disorder (Fig. 
3(a) I and when the disorder potential is stronger (Fig. 
3(b) I for fj, = 15. Co mpari son betwe en le ft plot of Fig. [T] 
(lower row) and Fig. 3(a) and Fig. |3(b)| shows that the 
strong localization of condensate density in co-ordinate 
space at /X = 15 is more due to the non-linearity (interac- 
tion) and not due to the disorder. 



20 r 



(a) 



The right plot of Fig. 3(a) md Fig. |3(b)| again plots 
the corresponding fc-space density distribution. The fea- 
tures agree with the Fig. [2] where a intermediate disorder 
strength is used with the showing of of density peaks at the 
corner of the first Brillouin zone apart of those in the cen- 
ter. As we have discussed the Brillouin zone density peaks 
is partly due to the Bragg refelction caused by numeri- 
cal discretization, but its height also depends on strength 
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(a) For fi — 15 the left figure (a) gives co-ordinate space den- 
sity in the absence of disorder for the condensate, obtained after 
the time of flight ui^t = r = 10. The right figure (b) gives the 
corresponding density distribution in the momentum space. 




(b) Similar plots as in Fig. 3(a) but for a stronger disorder than 
one used in Fig. \l\The left figure (a) gives the spatial density 
profile of the condensate. Plot (b) again gives the corresponding 
momentum space density. The stronger disorder potential has 
mean 0, standard deviation is 1.15 and the cutoff wavelength is 
Ac = 0.31. All other parameters are same as Fig. 



3(a) 



Fig. 3. Time evolved density of the propagating condensate 
with and without disorder 



of non-linearity, disorder as well and changes for the slow 
condensate. 

In Fig. [4] we have plotted the time evolution of the 
density profile for the entire time interval starting from 
ujzt = to ujzt = 10. From the actual calculation in the 
spatial grid we have chosen points at regular interval since 
plotting all the points make the size of the file too large to 
be represented through such color plots. But we have pre- 
sented the data in a way such that this omission of these 
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Fig. 4. color onhne:The figures in column (a) and (b) gives the 
full space-time profile of the density distribution for the same 
disorder potential used in Fig. [1] The x and y axes are same 
for all the plots. The chemical potentials has been indicated 
on the right side of each row while the right column again the 
initial density profile contains the slow k components satisfying 
kLTF{t = 0) < 1. 



intermediate data points does not affect the overall den- 
sity pattern. The upper plot in the left column for ^ — 10 
( weaker non-linearity) shows that the cloud localizes, ex- 
pands and then again localizes and does not spread mono- 
tonically. This is one of the new findings in this paper. 
However we have checked that the edge part of the cloud 
shows a decaying envelope all the time. At increased non- 
linearity for /i = 15 there is no such oscillation between 
spreading and localization and the cloud get strongly lo- 
calized in the lower plot of the left column. However the 
corresponding plot in right column where we remove the 
faster components from initially such stronger localization 
at /i = 15 is abscent and we get a non-monotonic behav- 
ior similar to the one for fi = 10. Thus we can argue that 
the lower Fourier components or equivalently the large 
wavelength components of the initial density is primarily 
responsible for this non-monotonic spreading of the prop- 
agating condensate. 

To summarize, in the preceding discussion we have 
studied the relative impact of non-linearity and disorder 
on the slow modes of initial density profile of the conden- 
sate by separating them from the faster ( almost ballis- 
tic) moving ones. We found that non-linearity has much 
stronger effect on the time evolution as compared to the 
effect of disorder. The later can be enhanced either by 



lowering /i or by removing higher k components from the 
initial condensate density. Now we extend this analysis for 
the time evolution of the phase and coherence properties 
of the condensate. 



3.2 Phase modulation, uncertainty and coherence 
properties 

It is known that in addition to the above described time 
dependent density, during time evolution the condensate 
attains a dynamical phase [U. In Fig. Q we plot this 
phase of the condensate as a function of position 
and the dimensionless time LUzt- The upper row of Fig. 
|5(a) compares the space time profile of the condensate 
phase for a given disorder but varying nonlinearity. As /i 
is increased from 6 to 20 to enhance interaction or non- 
linearity, the phase fluctuates on a shorter time and length 
scale. The details of the phase fluctuation and its degree 
of non-uniformity in space and time is however dependent 
on the particular realization of the disorder. For example, 
in the upper plot of column (b), at later time the phase 
fluctuation is more on the right side of the origin for the 
given realization of the the disorder. It could have hap- 
pened on the left side as well or at somewhere else for 
another realization of the disorder potential. The same 
plot ( first row, second column) also shows that the phase 
fluctuation is initially concentrated in the center and grad- 
ually spreads away. This happens as the initial profile of 
the static condensate spread away over a larger region in 
space with time and happens for any realization of the 
disorder potential. 



The lower row of the Fig. 5(a) gives the correspond- 
ing plots for the slow condensate. The removal of the 
faster Fourier components effects the interaction domi- 
nated phase fluctuation ( column (b) of the plots) differ- 
ently as compared to the disorder dominated phase fluc- 
tuations (column (a) of the plots). The interaction dom- 
inated phase fluctuation gets reduced in this space-time 
plot for the slow condensate where as in the left column 
disorder dominated (lower fi) phase fluctuation gets en- 
hanced for the slow condensate. However, this change in 
the fluctuation does not happen uniformly in the space- 
time. 

In Fig. 5(b) for = 6, we plot the space time profile of 



the condensate phase in the complete absence of the dis- 
order ( left column) and also in the presence of a stro nger 
disorder as compared to the one used in Fig. 5(a) The 



Fig. |5(b) as well as the upper figure of the left column 



Fig. 5(a) all of which correspond to /x = 6, shows that the 
introduction of the disorder potential and its variation has 
not any significant effect on the phase fluctuation. Thus 
increased non linearity has a stronger effect on phase fluc- 
tuation than increased disorder strength. 

Qualitatively, spatial phase fluctuations can be under- 
stood in the following framework. We know after time of 
order cozt ~ 1, most of the interaction energy gets con- 
verted into the kinetic energy. We can write the conden- 

i<t>u 



sate wavefunction as ip^ — ^J'pZi^ , then to a good ap- 
proximation ( if the disorder potential is much less than 
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the chemical potential) 



dt 



1 . dp^ 
Ap^ dx 



Since the typical density is lower than that of the static 
condensate, the most part of this kinetic energy is con- 
tributed by the phase fluctuation. Thus we can relate the 
high phase contrast region with the low density regions of 
the condensate. We can explain some of the regions with 
higher or lower phase contrast in this way. However this 
simple argument does not explain fully all the interesting 
features of the phase evolution under the combined effect 
of non-linearity and disorder. 

In the next figure Fig. [6] we plot the uncertainty in 
the position and the momentum in the cloud. They are 
defined in the following way. 



Ax = V a;2 - x^, 



J dxx"-p^{x) 
J dxp^{x) 



(4) 



where we have characterized the spatial distribution of 
the cloud by its moments. Similarly the uncertainty in 
momentum (fc) can we written as 



Ak= y B -k , A:" = 



Jdkk^pUk) 
J dkp^(k) 



(5) 



where we have characterized the fc-space distribution of 
the cloud by its moments, Axsl a-nd AksL £^re the cor- 
responding uncertainty in x and k when the initial con- 
densate is depleted of fast Fourier components such that 
kLxpit = 0) < 1. One of the purposes of these plots 
is to find out if there is any complementary relation be- 
tween Ax and Ak dictated by the uncertainty principle. 
For example, how much this behavior deviates from a co- 
herent state? Also the time dependence of Ax or Axsl 
tells us the nature of the condensate motion. For exam- 
ple, if Ax{t) ^ \ft the motion is diffusive, if Ax(t) ^ t 
the motion is ballistic. The saturation of the Ax{t) on the 
other hand implies inhibition of transport [15) . Thus calcu- 
lation of Ax{t) will provide a direct measure of the impact 
of non-linearity and disorder on the transport properties 
of the condensate. 

The plots (a) and (b) in Fig. 



6(a) are for the same 



disorder potential that has been used in Fig. [TJ Comparing 
the plots for Ax and Axsl in these two figures we find that 
the removal of the higher fourier components leads a more 
rapid temporal change in uncertainty. This can be clearly 
seen for a range of non-linearity ( Fig. 6(a) plots (a) and 
(h) , p ~ 10, 15). These results agrees with the observation 
in Fig. [1] where a strong localization in the special profile 
in the density for p — 15 is resisted and the condensate 
is allowed to spread over time in when the higher fourier 
components are removed from the initial profile of the 
condensate. Also the oscillator y time evolutions of Ax in 
Fig. |6(a) (a) and Axsl in Fig- 6(a) plots (a) and (b) are 
consistent with the time evolution of the density profile 
given in Fig. 01 Additionally these figures show that in 
course of such oscillatory behavior when the condensate 



spreads, the motion is almost ballistic since it is almost 
linear in time. 



A com parison between the plots Fig. 6(a) plot (b) and 
Fig. 6(b) plots (a) and (b) shows that the increase in the 
strength of the disorder does not affect the uncertainty 
very strongly. Also for p = 15 (higher non-linearity). Ax 
shows saturation after some time of fiight. We have checked 
that for even stronger non-linearity (such as p = 20) this 
saturation takes place in less times as compared to the 
cases of weaker non-linearity. 

plot 



Comparing plots of Ak and AksL in Fig. 6 (a 



(a) and plot (b) we find that as p is being increased, 
the uncertainty in k shows more fluctuations However we 
have checked that at even higher non-linearity ( for ex- 
ample at /i = 20) this fluctuation subsides. Comparing 
plots of Ak and AksL in Fig. 6(a) plot (b) and in Fig. 
6(b) columns (a) and (b) we also find that uncertainty 



in k fluctuates more with the introduction of the disorder 
potential. These fluctuations however subsides for the cor- 
responding slow condensate. Thus we may conclude that 
the removal of higher fourier components increases the 
uncertainty in position of the condensate and reduces the 
uncertainty in the momentum of the condensate. 

In the next two figures, namely Figs. [7] and [8] we plot 
the coherence in the evolving condensate as a function of 
space-time. To this purpose following |25| we have defined 
the function 



X{d,t) 



dxip^ix + d,t)^^{x,t)\^ 



(6) 



If x{d, t) — ^ there is full phase coherence and if x{d, t) = 
0, the phase coherence is said to be lost. For a given real- 
ization of the disorder potential and for a given value of d, 
the coherence function is obtained by summing over the 
product integrand in ([6]) for all pair of points on the dis- 
crete grid which are separated by this distance d at each 
time. If we assume that the ergodic limit is reached this is 
equivalent to the disorder averaging. We have also numer- 
ically verified that features of these plots are common to 
various realizations of the disorder potential. As the plot 
shows at a given time coherence changes as a function of 
the distance between two points. And for a given distance 
between two points, the coherence changes as a function 
of time. We shall discuss below the impact of the interac- 
tion as well as disorder on this spatio-temporal behavior 
of this coherence function given in ([6]) 

In Fig. 7(a) we plot such coherence as a function of 



space-time for different strengths of nonlinearity by chang- 
ing// but keeping the disorder potential same. The Fig. 
7(b) gives corresponding plots for the slow condensates. 



The Fig. 7(a) shows that for a given separation d between 
two points on the condensate profile, the coherence de- 
creases with time initially and then gets revived again. 
Such loss and revival of phase coherence has also been 
reported in numerical simulation of the BEC in a mov- 
ing optical lattice [25]. In Fig. 7(a) as we change from 
= 10 — > 15, with increasing non-linearity, the coherence 
decreases. Also the oscillation between loss and revival be- 
comes faster ( very thin stripes of alternating color in the 
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plot (b) of Fig. 7(a)). At a even higher non linearity at 
/z = 20, the coherence decays to a very low value at the 
very initial stage of the time of flight and the loss and 
revival behavior also disappears. In Fig. |7(b)[ for the cor- 
responding slow condensate both for /i = 15 and /i = 20 



the coherence is higher as comapred to those in Fig. 7(a) 
For /i = 15 the loss and revival feature now happens at 
a slower temporal rate and over a larger region of space. 
However in Fig. |7(b)| for /i = 20 the coherence decays 
to a very low value only after r 5 as compared to Fig. 
7 (a) I and there is no loss of revival behavior. These features 
imply that at a relative high value of chemical potential 
(/X = 20) the phase coherence decays as a function of space 
(-r^) and time and almost vanishes both in the case of 
large separation between two points at all times and also 
after long time over the entire length of the condensate. 

We contrast this situation with the cases with stronger 
disorder as well as the case without any disorder in the 
next Fig. [51 We consider a weakly non-linear case in Fig 
8(a) fo /i 
for 



6 and strongly non- linear case in Fig. 8(b) 
20. Comparison between left and right column 



of both figures shows that the coherence properties are 
almost unaffected by disorder irrespective of the strength 
of th e non-linearity. However a comparison between Fig. 
8(a) and Fig. 8(b) shows that that the coherence function 



indeed gets affected by the increase in non-linearity and 
the time evolved condensate becomes less coherent with 
increasing non-linearity. 

Within our purely numerical framework a complete 
theoretical understanding of these two intriguing features 
of coherence function, namely the loss and revival of this 
coherence for moderate non linearity ( such as the case for 
/Li = 10 or for slow condensate at /i = 15), and its inde- 
pendence on the disorder strength, is very difficult. Thus 
we confine ourselves in making some relevant comments. 
The integrand of the coherence function is dependent on 
both the density and phase modulation as — y/pZe^'^'^ ■ 
Using this the loss and revival feature can be related to 
to the successive spreading and contraction of the conden- 
sate density depicted in Fig. H) Generally we find that the 
coherence decreases for a spreading condensate the coher- 
ence increases for a contracting condensate. Now let us 
recall our previous discussion on phase fluctuation. Since 
the spreading of the condensate reduces its mean den- 
sity, the phase fluctuation increases. Increased fluctuation 
in phase leads to the loss of coherence. Similarly contrac- 
tion increases the mean density and thus suppresses the 
phase fluctuations w hich in turn le ads to the gain in coher- 

and Fig. [3(b)] 



ence. Also from Fig. 3(a) 



when the effect of disorder is enhanced in an otherwise 
non-linearity dominated bulk region by removing higher 
Fourier components from the initial profile of the conden- 
sate, it just modulates the condensate density according 
to the modulation of the disorder potential. However the 
envelope of the density profile does not carry the signa- 
ture of the disorder except at the edge where the density 
is very low. Our results suggests that because of this disor- 
der induced modulation the integrand in relation © will 
be changing very fast, the sum effect of this disorder in- 



duced density modulation will be almost zero. However 
this does not rule out the possibility that the higher order 
correlation functions involving tpw{x,t) will also remain 
insensitive to disorder. 

Thus our plots reveal that the spatio-temporal behav- 
ior of the coherence has more to do with non-linearity 
rather than disorder. A scattering theory based approach 
to this problem of loss of superfluidity in presence of dis- 
order has been considered recently in rcf. 26J. 



3.3 A brief comparison with a two-dimensional 
condensate condensate 

In this section we compare very briefly the above dis- 
cussed one dimensional cases with comparable situations 
in two dimension. To this purpose, for comparable chemi- 
cal potentials, we have plotted the density profile of two- 
dimensional condensates in a disorder potential after prop- 
agating it for dimensionless time 5 in position as well as 
in the momentum space. Also the chemical potential is in 
the unit of huj_i_ where lu± is the confinement frequency in 
the x-y plane. The results are plotted in Figs. [H and [TOl 
The initial profile of the condensate is again taken as an 
Thomas-Fermi condensate with parabolic density profile. 
Typically the condensate spreads more in two-dimension 
which agrees with the general assumption that in higher 
dimension the effect of disorder is less. The another in- 
teresting feature is that in the momentum space density 
profile we again see a localization of density in the Bril- 
louin zone boundary. Also after the filtering out of the 
faster components from the initial condensate profile here 
increases the spread of the condensate, which one can 
see particularly in the Fig. [TUlby comparing the two figures 
in the first column. A more detailed study of the conden- 
sate profile in two dimension in presence of non-linearity 
and disorder will be presented in a later submission. 



4 Summary of the effect of non-linearity and 
disorder on the condensate 

The main results of the preceeding discussion can be pre- 
sented as follows 

— We have tried to differentiate between the effet of dis- 
order and nonlinearity on the time evolution of the 
condensate in a regime where both of them are impor- 
tant. We found that non-linearity always overwhelms 
the effect of disorder 

— However the effect of non linearity can be reduced ei- 
ther by lowering the chemical potential or by removing 
faster Fourier components from the initial condensate 
profile 

— If the non linearity is not too strong and the effect of 
disorder is appreciable then the spread of the conden- 
sate goes to a non monotonic behavior by expanding 
and contracting during propagation recursively. This 
behavior is also accompanied by strong fluctuation in 
phase, higher uncertainty in momentum, and recursive 
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gain and loss in the coherence. With the increase of in 
the non hncarity this behavior is replaced by a strong 
localization of the condensate density in space. 

Thus altogether our numerical results show a very in- 
triguing non linear dynamics of the density and phase of 
the condensate when the effect of non-linearity and disor- 
der are present simultaneously . 

However all the above results are presented on the ba- 
sis of numerical work. The conclusions from such calcu- 
lations are always limited due to finite size effect, trun- 
cation of the evolution operator, numerical discretization 
etc. Nevertheless, they bring forth some very interesting is- 
sues and provide some insight for future work. Particularly 
the indication of the non- linear dynamics is intriguing and 
demands analytical treatment for a complete understand- 
ing. As we have pointed out in the beginning that in the 
successful experimental demonstration of Anderson local- 
ization [T^ for an interacting BEC, the density profile of 
the edge of the condensate is analyzed where effect of the 
non-linearity is negligible. Here we have mostly focused 
on the bulk region where the non-linearity is important 
to understand its effect on disorder induced quantum lo- 
calization. Our results will hopefully help to design exper- 
iments to probe the effect of non-linearity on Anderson 
Localization of matter waves. 
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(a) coior online: The uncertainty m position and momentum (k ) 
space as a function of time for two different values of fi for the 
disordered condensate. The corresponding /i has been indicated 
in each plot. Axsl and AksL are the corresponding uncertainty 
in position and momentum space for the slow condensate. The 
disorder potential is same as one used in Fig. Q) 



(a) color online: Phase (color axis) of the disordered condensate 
IS plotted as a function of space ( ,x-axis) and time (u^t, 
y-axis ) in presence of the disordered potential for fi — 6, (left 
column) and for /i = 20 (right column), x and y axes are same in 
all the plots. The lower row again gives such phase plot for the 
corresponding "slow condensate" for which kLTpit = 0) < 1. 
The disordered potential used is same as one in Fig. Q] 



) [1:15, no disorder 



(b) [1:15, strong disorder 
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(b) color online: Similar space-time plot of the condensate phase 
for /i = 6. Plot (a) is in the absence of disorder and the plot (b) 
is in the presence of a stronger disorder than one used in Fig. 



(b) color online: Similar uncertainty plot as in the Fig. 6(a) 
fi = 15. Left plot is for the case without any disorder and the 
right plot is for the case with a stronger disorder. The stronger 
disorder potential is same as the one used in Fig. The marker 
used for different plots are same as one used in Fig 



6(a) 



5(a) and same as one used in Fig.^ 



Fig. 6. Time evolution of the uncertainty in position and mo- 
mentum of the condensate 



Fig. 5. Time evolution of the phase of the propagating con- 
densate with and without disorder 
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(a) H^IO 



(b) H^15 (c) yu=20 

30 1 




(a) (1=6, no disorder 



(b) ^=6, stronger disorder 




5 10 

(a) color online: Coherence function x{d,t) (see Eq.\S^ is plotted 
along the color axis. The x axis is dimensionless time uj^t and 
the y axis is the parameter d in the scale of the Ltf of the 
corresponding condensate, fi is indicated in each upper plot. The 
disorder potential is same as one used in Fig. Ql Lower row gives 
plots for the corresponding slow condensates. 
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(a) color online: Plot of x{d,t) (color axis) as a function of the 
dimensionless time uj^t (x-axis) and (y-axis) for fj, — 6. The 
upper row is for the condensate. Left plot (a) corresponds to the 
case without any disorder and the right plot (h) corresponds to 
the case of stronger disorder potential previously used in Fig. 



a) [1=20, no disorder (b)n=20, strong disorder 



(b) color online: Similar plot of x{d,t) as in the Fig. \7(^ for 
the corresponding slow condensates 

Fig. 7. Spatio temporal evolution of the coherence of the con- 
densate for different nonlinearities (fi) 
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(b) The upper row is a similar plot of x{d,t) as in Fig. 8(a) 



hut for fj, — 20. In the lower we have plotted the corresponding 
behavior of the slow condensate 



Fig. 8. Spatio temporal behavior of coherence in the absence 
of disorder and in a stronger disorder 
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Fig. 9. color online: The density plot in co-ordinate and mo- 
mentum space with ^ = 10 after dimensionless time = 5. The 
disordered potential used here have mean 1.99 and standard 
deviation 0.23. (a) and (b) gives the density in x — y and kx, ky 
space for the full condensate, (b) and (d) gives the correspond- 
ing density plots after filtering out the higher momentum com- 
ponents at t = according to |A;|Ltf < 1- Here Ltf is actually 
Thomas Fermi radius in two dimension and |A:| — ^fc^ + ky. 
The length is measurd in the unit of Ltf and k is measured 
in the unit of L^p 




Fig. 10. color online: The density plot in co-ordinate and mo- 
mentum space with ^ = 30 again after dimensionless time = 5. 
The disorder potential is same as one used in previous figure, 
(a) and (b) gives the density in x ~ y and kx, ky space for the 
full condensate, (b) and (d) gives the corresponding density 
plots after filtering out the higher momentum components at 
t — according to \k\LTF < 1. Here Ltf is actually Thomas 
Fermi radius in two dimension and |A:| = ^/k^ + ky. The length 
is measured in the unit of Ltf and k is measured in the unit 
of L^p 
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